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ABSTRACT 

Feasibility of the geometric test as a probe of the cosmic equation of state of the 
dark energy is discussed assuming the future 2dF QSO sample. We examine sensitivity 
of the QSO two-point correlation functions, which are theoretically computed 
incorporating the light-cone effect and the redshift distortions, as well as the nonlinear 
effect, to a bias model whose evolution is phenomenologically parameterized. It is 
shown that the correlation functions are sensitive on a mean amplitude of the bias and 
not to the speed of the redshift evolution. We will also demonstrate that an optimistic 
geometric test could suffer from confusion that a signal from the cosmological model 
can be confused with that from a stochastic character of the bias. 

Subject headings: cosmology: theory - dark matter - large-scale structure of universe 
- quasars: general 

1. INTRODUCTION 

Recent observations of the cosmic microwave background anisotropics and the distant type la 
supernovae favor a spatially flat universe whose expansion is currently accelerating (de Bernardis 
et al. 2000; Lange et al. 2000; Perlmutter et al. 1999a; Riess et al. 1998). Motivated by this fact, 
variants of the cold dark matter (CDM) model with the cosmological constant (ACDM) have been 
widely studied. In particular, quintessence is proposed as a plausible model (e.g., Caldwell et 
al. 1998; Zlatev et al. 1998). An attractive feature of the quintessence is that certain quintessence 
models (tracker models) naturally explain the 'coincidence problem', the near coincidence of 
the density of matter and the dark energy component at present (Zlatev et al. 1998; Steinhardt 
et al. 1999, and references therein). Observational constraints have been investigated on the 
quintessential cold dark matter (QCDM) model (e.g., Efstathiou 1999; Perlmutter et al. 1999b; 
Wang et al. 2000; Newman &: Davis 2000), including attempts for reconstructing the cosmic 
equation of state or the quintessential potential (Starobinsky 1998; Saihi et al. 2000; Nakamura & 
Chiba 1999; Chiba & Nakamura 2000). 

In the meanwhile, recent progress on the 2dF QSO redshift (2QZ) survey has been reported by 
Shanks et al. (2000). An interesting scientific aim of the 2QZ survey is to obtain new constraints 
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on the cosmological constant from the geometric test, which was originally proposed by Alcock 
& Pancynski (1979), with the QSO clustering statistics. Several authors have proposed possible 
cosmological tests using the clustering of high-redshift objects (Ryden 1995; Ballinger et al. 1996; 
Matsubara & Suto 1996; Popowski et al. 1998; Nair 1999). In a series of work (Suto et al. 2000; 
Yamamoto et al. 2000, hereafter Paper I; references therein), a useful theoretical formula has been 
developed for predicting the correlation functions of high-z objects, incorporating the light-cone 
effect and the redshift-space distortions, as well as the geometric distortion, simaltaneously. The 
theoretical formula is useful because it is expressed in a semi-analytic form, which enable us to 
compute the two-point correlation functions corresponding to a survey sample based on specific 
cosmological models without huge numerical simulations. Therefore it will be worth to examine 
the feasibility of the geometric test as a probe of the cosmic equation of state of the dark energy 
with assuming a realistic QSO clustering statistics. 

Unfortunately, however, the QSO clustering bias is not well-understood at present. This 
ambiguity of the clustering bias should be crucial for the geometric test because the evolution 
of bias affects the predicted correlation functions. In Paper I, behavior of the QSO correlation 

functions is partially examined. On the other hand the stochastic bias has been discussed as a 
possible character of the galaxy biasing by several authors, (e.g., Dekel & Lahav 1999; Tegmark & 
Peebles 1998; Taruya et al. 1998). If the QSO clustering bias possesses the stochastic character, 
it affects the redshift-space correlation functions (Pen 1998). In the present paper we consider 
the feasibility of the geometric test as a probe of the cosmic equation of state of the dark energy 
component, including the possible stochasticity of the QSO clustering bias. This paper is organized 
as follows. In section 2, we briefly review the theoretical formula for the two-point statistics. In 
section 3, sensitivity of the correlation functions to the evolution of a bias model is investigated. 
Section 4 is devoted to summary and conclusions. Throughout this paper we use the unit in which 
the light velocity c equals 1. 



2. Contents for Theoretical Prediction 

We restrict ourselves to a spatially fiat FRW universe, and follow the quintessential 
cosmological model consisting a scalar field slowly rolling down its effective potential. The effective 
equation of state of the dark energy, wq = Pq/pq, can be a function of rcdshift in general case, 
however, we will consider the QCDM model with a constant equation of state for simplicity. In 
this case the dark energy density evolves pQ oc a{z)^^^^^'^^\ where a{z) is the scale factor. Then 
the relation between the comoving distance and the rcdshift is 

( sj_ r ^ ^ 

HoJo[n^{i + z')^ + nQ{i + z'f(^+^Q)Y/^' 

where Hq = 100/ikm/s/Mpc is the Hubble constant, fi^ and ^q{= 1 — ^m) denote the density 
parameters of the matter component and the dark energy, respectively, at present. 
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Wang & Steinhardt (1998) have given a useful approximate formula for the hnear growth 
index, which we adopt for computation of the /-factor defined by f{z) = d\nDi{z)/dlna{z), 
where Di{z) is the hnear growth rate. Simple fitting formulas for the linear transfer function and 
the nonlinear mass perturbation power spectrum are presented by Ma et al. (1999) for the QCDM 
model, which we adopt for modeling the CDM mass density perturbations. The formulas are 
applicable for the QCDM model with the constant equation of state — 1 < wq ,^ — 1/3 and up to 
2; ~ 4 with errors ^ 10 %. In addition we adopt the normalization by cluster abundance (Wang 
&; Steinhardt 1998; Wang et al. 2000). Throughout this paper we assume the Harrison-Zeldovich 
spectrum. 

In predicting clustering statistics of high-redshift objects in a redshift survey, several 
observational effects must be incorporated for careful comparison between theoretical predictions 
and observational results. A useful theoretical formula for the two-point statistics has been 
developed incorporating the redshift distortions due to peculiar motion of sources and the 
light-cone effect simaltaneously, as well as the geometric distortion (Suto et al. 2000; Paper I). 
According to the theoretical formula, the two-point correlation function is given by a Fourier 
transform of the power spectrum on a light cone Pi"^{k) 

1 

ii{R) = dkk^MkR)Pt''ik), (2) 

where P^^^(/c) is obtained by averaging the local power spectrum Pf'^^{k,z) over the redshift 

Lc.> ^ _ !dzW{z)Pr{hz) 
JdzWiz) 

with the weight factor 

/dNy/s^ds^ 



where dN/dz denotes the number count of the objects per unit redshift and per unit solid angle, 
and s = s{z) denotes the distance-redshift relation of the radial coordinate that we chose to plot 
a map of sources. In the present paper we adopt the distance-redshift relation of the Einstein de 
Sitter universe, i.e., s{z) = 2Hq^[1 — (1 -|- z)"^/^]. In equation (^ z-integration arises from the 
light-cone effect within the small-angle approximation. The power spectrum Pi^^{k^z) in (^ is 
given by 

2/ -I- 1 

pr{k,z) = ^ / d^iLiiii) 

cic^\ Jo 

xPqsoUii ^ — , |q±l ^ (5) 

V C|| Cx / 

where PQso{q\\,\(l±\, z) is the QSO power spectrum, (q^) is the wave number component 
parallel (perpendicular) to the line-of-sight direction in the real space, and Li{ii) is the Legendre 
polynomial. In equation (||), we denoted k = |k|, Cx = r{z)/s{z) and = dr{z) / ds{z) with the 
comoving distance in the real space r{z), equation (|l|). 
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We model the power spectrum of QSO distribution by introducing the bias factor b{z) 



PQso{q\\, |q. 




(6) 



where q = \/qf^+ and -Pmass(^i-z) is the CDM mass power spectrum. The terms in proportion 

to f{z) in (y) is traced back to the Hnear distortion (Kaiser) effect. To describe a cross correlation 
between the QSO distribution and the CDM mass distribution, we here introduced the cross 
correlation coefficient TZ{z). In the case of the deterministic bias, the cross correlation coefficient 
can be set 'JZ{z) = 1, however, in the case of the stochastic bias, the cross correlation coefficient is 
allowed to deviate from unity (e.g., Pen 1998). 

In (^), D[qi^ap{z)] is the damping factor due to the Finger-of-God effect, for which we adopt 
the exponential model for the distribution of the pairwise velocity dispersion ap. For ap we adopt 
an approximate formula whose validity is investigated by several authors (see Magira et al. 2000; 
Mo et al. 1997). The nonlinear effect is not a dominant effect on large length scales, however, it 
can not be neglected. Therefore we here take the nonlinear effects into account for definiteness. 

Recently, Boyle et al. (2000) have reported on the evolution of QSO optical luminosity 
function using their preliminary result of the 2QZ survey. The evolution of the luminosity function 
is compared with analytic fitting formulas. For theoretical predictions which correspond to the 
on-going 2QZ survey, we adopt their quasar luminosity function with best-fitted parameters of 
a power-low polynomial evolution of luminosity under the assumption of the Einstein de Sitter 
universe. Then the number count dN/dz brighter than the limiting magnitude B < 20.85 can be 
obtained by integrating the luminosity function. [| 



The QSO clustering bias is a challenging problem and a few authors have proposed theoretical 
models (Fang & Jing 1998; Martini &: Weinberg 2000; Haiman & Hui 2000). However, these models 
seem to be prototype models, hence we here adopt a model whose evolution is phenomenologically 
parameterized, for simplicity, as (c.f., Matarrese et al. 1997), 



where a, (3 and b{z^) are free parameters, /3 specifies the speed of redshift evolution, 6(2:*) is the 
amplitude at a mean redshift , and a corresponds to the amplitude of bias at z = in the limit 



For the K-correction we assume the quasar energy spectrum L^, oc v 



3. Sensitivity 




(7) 
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of 1- Throughout the present paper we fix a = 0.5, however, this does not change the results 
for ^ a ^ 1. We define the mean redshift by 

pmax (IZZW(Z) 

with the weight factor (^). Here we consider the sample in the range 0.3 < z < 2.2, in this case we 
have = 1.2. We note that is almost same as the mean amplitude of the bias defined in the 
similar way to @. 

For the geometric test, the ratio of the correlation functions (,2{R)/(,oiR) will be an important 
quantity to characterize the geometric distortion effect. Figure 1 shows contours of 5,2 /(,o on the 
(/? — b{z^)) plane with the separation R fixed as 20/i~^Mpc for various cosmological models {solid 
lines). Panels (a) and (b) show the cases of the ACDM model (Qm = 0.3, wq = — 1) and the 
QCDM model {^m = 0.3, wq = —1/2), respectively. In both panels we show the case of the 
deterministic bias. This figure shows that the value of is sensitive only to the parameter 

and not to the speed of evolution f3. 

Recently it is reported that the QSO correlation function is consistent with being 
^ = (r/ro)^^'^ with tq = 4 /i~^Mpc from the preliminary result of the 2QZ survey (Shanks et 
al. 2000). Previous result of the QSO surveys similar to the 2QZ survey reported rg = 6 /i~^Mpc 
(Groom & Shanks 1996). The dashed lines on the panels show the contours satisfying 
i?c = 4 /i~^Mpc and Rc = Q /i~^Mpc, where the characteristic correlation length Rc is defined 
by (,o{Rc) = 1. The observational constraint is not strict at present because of statistical errors. 
Then Figure 1 should be regarded as a demonstration. However, it is instructive. Taking the 
constraint 4 h~^Mpc < i?c ^ 6 /i~^Mpc into account, the ratio of the correlation function is 

- -2.0 1.4 for the ACDM model, and ^2/^0 - -1-4 ~ -1.0 for the QCDM model. 

We fixed as i? = 20/i"iMpc in Fi gure 1, however, similar behavior appears at other values of 
R. For example, in the case R = 30/i-^Mpc, 6/^0 - -3.5 ~ -2.3 for the ACDM model and 
(,2/^,0 — —2.3 ~ —1.8 for the QCDM model. The reason why (,2/^,0 depends on the cosmological 
model can be explained as follows. One reason is the difference of the linear growth rate Di(z), 
which affects the Kaiser factor. The other effect is the scaling effect due to the geometric 
distortion, which is described by the factors Cy and c±. The latter effect is dominant. It is clear 
from Figure 2 that the less negative value of wq {wq > — 1) decreases the coefficients C||(z) and 
c^{z). This causes the difference in depending on the cosmology. 

Now let us discuss a possible effect from stochastic character of the bias by considering the 
case TZ{z) 7^ 1 in (^). For simplicity we consider the case TZ{z) is constant. The panel (c) in Figure 
1 plots the contour of ^2/Co of the ACDM model same as the panel (a) but with lZ{z) = 0.6. 
As is clear from the panels (b) and (c), the both panel becomes almost same. That is, the ratio 
£,2/£,o takes the same values in the cases {wq = —1/2, TZ = 1) and {wq = —1, TZ = 0.6). Thus 
the signal from the cosmological model can be confused with that from the stochastic character of 
the bias. This fact suggests that an optimistic geometric test suffers from degeneracy between the 
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cosmological model parameter and the stochasticity parameter of the bias. 

The panel (d) in Figure 1 is same as the panel (b) but with = 0.35, which shows 
Om-dependence. By comparing (d) and (b) , ambiguity of the density parameter Qm might not be 
negligible too. 

To solve the degeneracy due to the stochastic character of the bias, ^4{R) might be useful. 
Figure 3 plots contours of ^4/^,0 for the various cosmological models, whose parameters are same 
as those in Figure 1. Similar feature to Figure 1 can be seen in Figure 3. However, Figure 3 shows 
that ^4/^0 is rather insensitive to the stochasticity TZ and the density parameter Qm, as is expected 
from the investigation of the linear stochastic biasing in redshift-space (Pen 1998). Therefore 
C4/C0 might be useful to break the degeneracy, if it could be measured precisely. However, the 
amplitude of the signal is rather small, which is of order of 10% compared with ^2/^0- 

4. Discussions 

In the present paper, we have examined the sensitivity of correlation functions in a QCDM 
cosmological model to the bias evolution, incorporating the various observational effect, i.e., the 
light-cone effect, the linear distortion effect, and the non-linear and the finger of God effects. 
Then the feasibility of the geometric test is discussed as a probe of the cosmic equation of state 
assuming the future 2dF QSO sample. The amplitude of the correlation functions is sensitive 
to the mean amplitude of the bias, and is rather insensitive to the speed of evolution due to 
the light-cone effect. We have found that the systematic difference appears in the ratio of the 
correlation functions depending on the effective cosmic equation of state, due to the geometric 
distortion effect. We have also shown that, if the QSO bias has a stochastic character, the signal 
from the cosmological model can be confused with that from the stochastic bias. Hence the simple 
geometric test with only ^2/^0 suffers from the degeneracy between the cosmological parameter 
and the bias parameter unless the stochasticity character is clarified. For example, the ACDM 
model with Tl{z) = 0.6 and the QCDM model with wq = —1/2 and 'R-{z) = 1 predict almost 
same value of at i? = 20/i~^Mpc. Therefore other cosmological information is required, e.g., 

the higher order multipole moment of the correlation function $.4:/(,o- However the signal from ^4 
seems to become noisy, more detailed investigations will be need for the viability. 

Finally it will be worth to discuss the robustness of our results for several assumptions adopted 
in the present paper. Another choice of the cosmological redshift-space s{z) alters the shape of 
the correlation function (Paper I), hence the predicted values of C2/C0 will be altered. However, 
the sensitivity on the cosmic equation of state will not be significantly altered, neither will be the 
feasibility of the geometric test. The bias model (equation [7]) seems to express general evolution 
of the bias, then our result is not sensitive to the bias model unless the scale-dependece of the bias 
is significant. Concerning the QCDM cosmological model, our investigation is restricted to the 
case Wq is a constant. And we used the fitting formula by Ma et al. (1999), which is applicable to 
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that case. In most cases, the quintessence equation of state changes slowly with time, however, we 
believe that predictions are well-approximated by treating wq as an averaged constant value (e.g., 
Wang et al. 1998). An open CDM model shows the similar result with the QCDM model in 
(Paper I), then our conclusion is based on the assumption of the spatially flat universe. 

We thank Y. Kojima, A. Taruya and T. Chiba for useful discussion and comment. This work 
is supported by the Inamori Foundation. 
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Fig. 1.— {Solid line) Contours of e,2{R)/^o{R) at R = on the (/? - plane 

for various cosmological models: (a) ACDM model ($7^ = 0.3, tag = —1.0); (b) QCDM model 
{Qm = 0.3, wq = —1/2); (c) same as the panel (a) but with the stochastic bias case with Tl{z) = 0.6; 
(d) QCDM model ($7^ = 0.35, 'u;q = —1/2). The contour levels of the solid line are indicated on 
the figure, and the same contour levels are adopted for each panel. In each panel we adopted 
h = 0.7, ^hh^ = 0.015, and a = 0.5 for the bias model. Except for the panel (c), the deterministic 
bias is considered. The left (right) dashed line shows the contour that satisfies the condition of the 
characteristic correlation length Rc = 4/i~^Mpc {Rc = 6/i~^Mpc). 
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Fig. 2. — Evolution of C||(z) and Cj_{z). The lines show the cases = 0.3 and wq = —1, —2/3, 
—1/2 and —1/3, from top to bottom, for both panels. 
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Fig. 3. — Contours of ^4(i?)/^o(-R) x 10 at i? = 20h ^Mpc on the (/3 — 6(2;*)) plane for various 
cosmological models. The model parameters and the meanings of lines are same as those of Figure 
1. 



